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D ■ Abstract 

\f~^ • We consider a class of plerions that have been detected in the TeV range, and 

. investigate the possibility that the emission is due to vr^ decay. From the TeV flux 

we derive what is the expected u flux at Earth and find that this is likely to exceed 
J> ' the detection threshold of the upcoming km^ neutrino detectors. 
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S ■ 1 Introduction 

Plerions are supernova remnants (SNRs) with a filled morphology. These rem- 
^ ' nants are characterized by a center-brightened nebula often seen in the radio 

■ and X-ray wavelenghts and believed to be powered by an embedded pulsar. 

Typified by the Crab Nebula, they have non thermal spectra at all wave- 
lengths. They have a flat power-law spectral index (a ~ — 0.3, 5*1^ ~ z/~") 
in the radio band and hard photon index (7 ~ 2,7 = a + 1) in the X-ray 
band. Out of ~ 220 Galactic SNRs only about 10% are classified as plerions 
(or Crab-like SNRs) [10]. 

Plerion spectra are usually well interpreted from the radio to the X-ray band 
as synchrotron emission of a population of relativistic pairs continuously sup- 
plied by the central pulsar. At higher energies the Inverse Compton Scattering 
of the same electrons and positrons off either an internal or external target 
radiation can play a role. However it is not clear whether this latter process 
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can be responsible for the emission recently observed from a few objects at 
TeV energies [16]. An alternative mechanism to produce TeV photons may be 
the decay of neutral pions produced through nuclear collisions of relativistic 
protons. 

In this paper we investigate the consequences of a possible hadronic origin of 
these TeV detections. Following the line of a recent work by Alvarez-Muhiz 
& Halzen [3], we compute the high energy neutrino flux at earth and find 
that the predicted fluxes may be detectable by large, km^ effective area, high 
energy neutrino telescopes, such as the planned south pole detector IceCube 
[12] or the Mediterranean sea detectors under construction (ANTARES, [4]; 
NESTOR, [14]) and planning (NEMO, [15]; see Ref. [11] for a recent review). 



2 TeV observations of plerions 

Four plerions have been so far detected at TeV energies, while upper limits ex- 
ist for a few others. The objects for which the detection is at a high confidence 
level ( ^ 4(t) are: the Crab Nebula [2], the Vela X SNR [20], the pulsar wind 
nebula around PSR1706-44 [13] and the radio nebula surrounding PSR1509- 
58 [16]. The VHE emission is unpulsed and therefore hkely to be associated 
to the pulsar wind nebula rather than to the pulsar magnetosphere. 

In Table 2 we report the list of pulsar wind bubbles which have been detected 
at TeV energies, supplied with the central pulsar luminosity and distance (Ref. 
[1] and references therein), and the observed TeV spectrum. 

Two of the objects in the table deserve some comment. First of all, it should be 
noticed that the association between pulsar B1706-44 and the remnant G343.1- 
2.3 is questionable as discussed by Giacani et al.[9], and in the following we 
refer to the radio nebula detected by Frail et al.[6] as the remnant associated 
to this pulsar. 

As to B1509-58, this pulsar is found in a very extended supernova remnant 
with a complex morphology. However a synchrotron nebula has been found 
with confidence surrounding the pulsar at X-ray frequencies [17,18,5], although 
no pulsar wind bubble has been detected at radio frequencies [8]. Moreover 
the spectral index at TeV energies has not been determined with confidence 
[16]. In the following we use a value of 2.5 in analogy with the spectra of the 
other objects and derive the normalization from the integrated photon flux 
measured by CANGAROO. 

As we mentioned in the introduction, a possible mechanism to interpret the 
TeV fluxes is the Inverse Compton Scattering (ICS), on the ambient photon 
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Table 1 

Pulsar wind bubbles detected at Te V energies. The name of the pulsar and its associ- 
ated remnant are reported in the first two columns. The pulsar bolometric luminosity 
and distance from Earth are in the third and fourth column respectively. In the last 
column the Te V spectrum as given in the references cited above is reported. 
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field, of the electrons responsible for the synchrotron emission at lower fre- 
quencies. The target radiation for ICS is given by the sum of two different 
contributions: the external photon field and the internal radiation produced 
by the source itself. A detailed discussion of the relative importance of the dif- 
ferent contributions, within the context of the objects considered here, can be 
found in Ref. [1]. The conclusion is that the target photon density is expected 
to be dominated by the cosmic microwave background radiation (CMB) at 2.7 
K, corresponding to an energy density Wph = 0.26 eV cm~^. 

When one considers the ICS and synchrotron emission in frequency ranges 
such that they are produced by the same electrons, the ratio of the syn- 
chrotron and ICS luminosities {Lgync and Lies, respectively) is equal to the 
ratio between the energy density of the magnetic and radiation fields. Since, 
as we mentioned above, the target radiation is mainly made of CMB photons, 
we can write, for the ratio between the ICS and synchrotron emission: 

-^ICS _ WCMB ^-^^ 
-^sync "^B 

where wcmb and wb indicate the CMB and magnetic field energy densities 
respectively. Eq. 1 can be used to estimate the magnetic field, once we specify 
the frequency bands in which Lies and Lgync have to be taken. 

The energy Eg of an electron producing TeV photons by upscattering the 
CMB is given by: 

E,[TeY] ^ 20 (;^) , (2) 
where E^ is the photon energy. Therefore the synchrotron photons emitted by 
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the same electrons will have a typical energy of: 

e,ync - 0.08 5_5 keV , (3) 

where is the nebular magnetic field in units of 10^^ Gauss. This value 
of the magnetic field strength is of order of that typically estimated for these 
nebulae assuming equipartition. 

If the TcV flux is due to ICS on the CMB, we then find for the magnetic field 
strength: 

/ L \^/2 

S/c = 3 X 10-^ — ^ Gauss . (4) 

Vi^TeV/ 

where Ltcv is the luminosity at TeV energies obtained integrating the TeV 
spectrum given in Table 2 above 1 TeV. 

In Table 2 we compare the value of Bjc with what can be found assuming 
equipartition B^q. For the flrst three remnants we have computed the equipar- 
tition magnetic field based on the radio data while in the case of the nebula 
surrounding PSR B1509-58 we have used the luminosity and size at X-ray 
frequencies, as given in Ref. [17]. 

Table 2 

Values of the parameters used to estimate B^q and Biq. The computed fields are 
reported in columns 6 and 8 respectively. In the second column we report the radio 
luminosity integrated over the frequency hand specified in the adjacent column. The 
spectral index and the extension of the nebula at radio frequencies are reported in 
columns 4 o-nd 5 respectively. The numbers we report for Crab and Vela are taken 
from the SNR Catalogue [10], while those for the nebulae associated to B1706-44 (^'^-d 
B1509-58 are from the references given in the text. Finally, the X-ray luminosities 
given in column 7 are taken from [1]. 
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It is apparent from the table that while for two of our four sources the ICS 
and the equipartition values of the field are very close, in the case of B1706-44 
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and Vela a noticeable discrepancy (a factor of 4 and 10 respectively) is found. 
For these two objects the possibility that the TeV emission is due to hadronic 
processes is particularly appealing. If this is actually the case, then we expect 
a neutrino flux from these sources that we compute in the next section. 



3 Neutrino events 



A way to disentangle electromagnetic and hadronic sources of high energy 
7-ray emission observationally is to look for the neutrino signals. Relativistic 
protons may produce TeV 7-rays either by photo-meson production or inelastic 
nuclear collisions. The relative importance of the two processes depends on the 
target density of radiation and matter in the source. The main difference, as 
far as their outcome is concerned, is in the fraction of energy that goes into 
charged pions compared to neutral ones. This translates in a different ratio 
between the total neutrino and photon energy flux. In the case of plerions 
the most likely process at work is p-p scattering, as can be readily seen by 
comparing the rates of photomeson production and p-p scattering estimated 
below. 

For photo-meson production the target for high energy protons is the plerion 
emission. The fractional energy loss rate of a proton with energy Ep (— F rUp (?) 
due to pion production results in (Ref. [19]): 
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Here we have treated the plerion as homogeneous and used the fact that the 
photon spectrum is a power law. F^, oc u'^. We have also made the approx- 
imation that the main contribution to pion production comes from photon 
energies e^, ?a epcak=0.3 GeV, where the p-7 cross section peaks due to the A 
resonance. The numerical values are obtained using: cXpeak = 5 x 10^^^ cm^, 
Cpeak =0.2, Ae=0.2 GeV, Up = epeak/(r h) and Pp ~ 1. Finally we have scaled 
the nebular radius and distance to the typical values of 1 pc and 1 kpc, re- 
spectively ( Rpc = Rat/pc and dkpc = d/kpc), and expressed the nebular 
synchrotron radiation flux in units of mJy. 

The energy loss-rate of a relativistic proton due to inelastic nuclear coUisions 



5 



can be estimated as 

.-1 , t Miv 3 

~ C t^o c C 1— Bs" ^0 c , (6) 



m. 



■p 47rR^ 



where nt is the target density, which we have expressed in terms of the nebular 
radius Rn and content of thermal material M^. Introducing in Eq. 6 the 
numerical values of the cross section for p-p scattering, (Tq = 5 x 10~^^cm^, 
and of the average fraction of energy lost by the proton, ( ~ 20%, we obtain: 

» 10-' ^ yr-' . m 



where M-^q = M^/Mq. From the comparison between Eq.7 and Eq.5, it is 
apparent that nuclear collisions are by far the most likely mechanism for pion 
production in plerions. 



From this follows that the relation between the neutrino and photon flux is: 



^1/ . Uf J-V^ 

IPmin Rmin 



where i^;™" (E'"^^^) is the minimum (maximum) energy of the p-p produced 
photons and E""" ^ (l/2)i?""° and E'J^"^^ ^ (1/2)^™^'' are the corresponding 
minimum and maximum neutrino energies. 

We estimate the neutrino flux in the energy range 1-100 TeV, which is the 
range in which the future km^ neutrino detectors will operate. Therefore we 
are concerned with photon energies above 2 TeV. 

In Table 3 we report the total number of neutrino events expected in a yr of 
operation of a km^ detector. This is obtained directly from: 

7V^ = 4AffxTx J j^{2E,) P,^{E,) dE, , (9) 

ITeV 



where P^^ = 1.3-10 ^ £^i/,TeV [7] is the detection probability for neutrinos with 
Ey ^ 1 TeV, T is the observation time and A^s is the effective area. 

The number of atmospheric neutrino events collected in a km^ detector during 
1 yr is of order 1. This is estimated assuming a background neutrino spectrum 
0t/,bkg ~ 10~'^-^i7TeV cm^^ s^^ sr^^ for E^ > 1 TeV, and a detector angular 
resolution of 0.3° hke NEMO [15]. 
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Table 3 

Predicted number of muon events, N^, in a krr? detector, using Eq.9 for the plerions 
considered in the paper. One year of integration time is assumed. 
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Our calculation suggests that the neutrino signal from three of the plerions 
considered here maybe identified above the atmospheric neutrino background 
by the next generation underwater (ice) neutrino telescopes. 



4 Discussion 

We have investigated the implications of a possible hadronic origin of the high 
energy emission from the plerions which have been detected by CANGAROO 
at TeV energies. The alternative explanation for this emission is based on ICS. 
Prom the comparison of the implied magnetic field strength under this latter 
assumption with the equipartition value, we have found that the ICS origin 
of TeV photons poses some difficulties for at least two of our four sources, 
namely Vela and the remnant surrounding B 1706-44. 

We have computed the implied neutrino flux under the assumption that the 
emission seen by CANCAROO is due to p-p scattering. We conclude that, 
for three of these objects, this hypothesis can be confirmed or disproven with 
confidence, given that the upcoming neutrino telescopes will be operated for 
a few years. 

The positive detection of neutrinos with Ei^ ^ TeV would not only provide 
the final answer to the origin of the detected VHE 7-rays but also confirm the 
hypothesis that protons can be accelerated upto at least tens of TeV in these 
objects. 
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